Anomalously fast kinetics of lipid monolayer buckling 
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We re-examine previous observations of folding kinetics of compressed lipid mono- 
layers in light of the accepted mechanical buckling mechanism recently proposed [L. 
Pocivavsek et al, Soft Matter, 2008, 4, 2019]. We set conservative limits on a) the 
energy released in the mechanical buckling process and b) the kinetic energy entailed 
by the observed folding motion. These limits imply a kinetic energy at least thirty 
times greater than the energy supplied by the buckling instability. The most likely 
source for the needed kinetic energy is self adhesion of the hydrophobic faces of the 
fold, not included in the accepted mechanism. 
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I. INTRODUCTION 



When a flat thin sheet of material is subjected to increasing pressure, it eventually buckles, 
crumples, or cracks. By buckling, the sheet releases stress over wavelengths much larger 
than its thickness. During the past years there has been increasing interest in phenomena, 
such as crumpling, where the deformation goes from an initial uniform state to a localized 
region occupying an arbitrarily small fraction of the sample plj. Monolayers of surface- 
active molecules (surfactants), adsorbed on a liquid, are found in many systems containing 
water-air or water-oil interfaces [2]. Such monolayers exhibit a rich variety of collapsed 
structures under lateral compression, most of which occur on an intermediate scale between 
the macroscopic one and the molecular thickness of the layer. Fluid monolayers collapse 
into disks, tubes, or pearls-on-string structures, depending on spontaneous curvature and 
charge of the lipid monolayer [3j. Like elastic sheets, many solid-like monolayers crack 
[U E] or buckle [6] under pressure, yet other solid-like monolayers fail by abrupt buckling 
into straight, micron-wide folds [7HT5] (See Fig. [l] for an example). In addition, liquid-like 
monolayers may form micron-scale vesicular objects of various shapes |71 |16l [17] or giant 
convoluted folds [U [TU] . This type of folding is believed to be driven by the interfacial energy 
gained from the contact across the two sides of the fold [TU] . Thus, the failure of surfactant 
monolayers under lateral pressure displays distinctive mechanical behaviors, which crucially 
depend on the in-plane rigidity [3 [HI [T9| and viscoelasticity [lO]. Lipid monolayers and 
bilayers are the material of choice for spatial partitioning in living matter, such as cells and 
compartments within them. These partitions are often observed to fold and wrinkle under 
stress. In particular, lipid monolayers that model the expanding and contracting sacs in an 
animal's lung exhibit the abrupt folding signature described in Refs. [2UH22]- 

This abrupt buckling is particularly puzzling. The abrupt relaxation motions or "jerks" 
have been extensively characterized for a particular system [13] , but their distinctive length 
and time scales have not been explained. Why do the jerks occur over a time scale of about 
0.1 seconds — far from any molecular or apparatus time-scales? Why does the motion stop 
suddenly, before the driving stress is relaxed? Why do the jerks show a robust characteristic 
displacement of a few microns? 

An intriguing hypothesis to explain this micron length scale was recently proposed by 
Pocivavsek et al [H [H] . These authors recalled that any thin sheet under compression on a 
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FIG. 1. The bright hne is a fold created in a compressed mixture of DPPG:POPG lipids. Dark 
domains are the condensed phase surrounded by a liquid phase in purple. Bar indicates 50 fim. 
Figure taken from Ref. |13j . 



liquid substrate buckles at a well-defined wavelength A. For lipid monolayers the expected 
value of A lies in the micron range. Thus it is of the same order as the characteristic 
displacement of the jerk relaxations. These authors also noted that the incipient wrinkles 
at wavelength A are unstable against folding, in which the excess wrinkled material from 
throughout the sample is concentrated into a single loop or fold. 

This paper aims to account for the dynamical features of the jerk motion. In Section [ITl we 
argue that the monolayer may be viewed as a broad, thin slab that translates almost rigidly 



during the jerk motion. Then in Section III we describe how such a slab should respond 



to horizontal forcing, accounting for progressive viscous entrainment of the fluid subphase. 
Using this result, combined with experimental observations, we set a lower bound on the 
kinetic energy of the jerks. We then survey the possible forces that might give rise to this 
kinetic energy, notably the energy released by folding in the mechanism of Ref. [11]. Even 
upper limits of this folding energy are far less than this observed kinetic energy. In Section 



IV we discuss effects that might account for this kinetic energy. 
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II. SYSTEM 



For definiteness we focus our study on the jerks analyzed in Ref. [13] and re-examined 
in Ref. pH]. We begin by reviewing the parameters of this system. The monolayer was 
a 7:3 mixture of dipalmitoylphosphatidylcholine (DPPC) and dioleoylphosphatidylglycerol 
(POPG) spread on pure water at room temperature in a 15-cm-long Langmuir trough and 
viewed in a microscope in a 100-micron-wide field of view at conventional video frame rates. 
It is compressed at 0.1 mm/sec to a nominal pressure of about 70 mN/meter before viewing. 
The measured pressures are consistent with the slight overpressure needed for folding [HI |23] . 
In these conditions this material microphase separates into a biphasic foam-like pattern 
of compact patches separated by narrow strips of different composition, visualized by a 
dilute fluorescent additive. The patches are 15-25 microns in size. The monolayer behaves 
mechanically as a solid, not a fluid; the jerks move in the direction of the Langmuir barrier 
motion. Rheological measurements in similar systems showed stress relaxation times of 
order 10-10^ sec. [TOl [2l] . The jerks vary statistically in their net displacement A and their 
duration t. No statistical correlations between jerks were observed. The displacements A 
vary from a minimum of about 1.2 microns to several times larger, with an average of 2 
microns [13j. The durations t vary from a minimum of about 0.09 sec to a few times longer, 
where the average is 0.12 sec. It has been suggested p^j that the larger jerks are cascades 
of elementary jerks, with this average A and t. In the estimates below, we use these average 
A and t values. 

We may simplify our description of the monolayer using three further features. First, the 
monolayer may be assumed to translate rigidly outside the folding region over distances of 
several cm. It thus entrains substrate fluid over these distances. In principle compressibility 
could invalidate this assumption. Compressibility implies that an initial imbalance of mem- 



brane stress produces a compressional wave whose speed is given by, c = yY/p^, where Y is 
the two-dimensional uniaxial compression modulus and ps is the mass density per unit area. 
If forces are applied on a timescale t, the resulting compression or expansion is confined to 
distances L < ct. Conversely, the compression is negligible and the body moves rigidly if 
its size L is smaller than about ct. Thus to show that our monolayer translates rigidly, we 
must establish that c is sufficiently large. 

For our case the density ps is the density of the material to be accelerated in propagating 
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the wave. Since the monolayer entrains substrate fluid as it moves, the density must take 
account of this subphase. We may flnd a lower limit on this c by using a lower limit to the 
modulus Y and an upper limit for the surface density ps- Accordingly we estimate Y by 
neglecting any compressional effects of the subphase. Measurements of the monolayer com- 
pression modulus range from 0.1 N/m to several N/m [18l [25], [26] . We use the conservative 
estimate of y > 0.1 N/m below. As for the surface density ps we flnd an upper limit by 
including all the water that might be entrained. As discussed in Section |III[ an upper limit 
on the entrained density ps is given by ps < 2p\/iH, where p {10^ kg/m^) is the density of 
the water subphase, and u (10~^ m^/sec) is its kinematic viscosity. Thus in the jerk time 
of 0.12 sec, a local compression can propagate a distance c t > 0.05 m, i.e. even by our 
minimal assumption, information from the folding region has reached about half the sample 
(5 cm) within the jerk time. Thus we expect only minor effects from compressibility of the 
sheet. In the calculations to follow we will assume L = 5 cm is the size of the moving sheet. 

A second feature of the monolayer is that its compressibility has negligible effect on the 
energetics of folding. As explained below, the over-pressure p released in the wrinkle-to- 
fold transition [H] in this system must be smaller than 6 x 10^^ mN/m. Given the large 
lower-bound modulus Y above, the compressive displacement Acom in a system of length L is 
smaller than Lp/Y < 0.03 micron. It is thus much smaller than the observed displacements 
A. Likewise, the compressive energy released is a small fraction (pAcom / pA) of the folding 
energy. 

A third feature of these monolayers is that the longitudinal propagation of the fold tips 
is much faster than the transverse folding and its accompanying jerk. Within the temporal 
resolution of Ref. (~ 0.03 s) the fold traverses the fleld of view (~ 150 pm) instanta- 
neously; this sets a lower bound of ~ 5 mm/s for its speed. In the system of Ref. [TT] two 
types of longitudinal folds were observed, the slower of which propagated at ~ 10 mm/s. 
Comparing these values with the characteristic translation velocity, 10 pm/s [13j, we see 
that the requirement is safely fulflUed. These features allow us to consider a simplifled two- 
dimensional problem of a thin elastic sheet moving over a semi-inflnite viscous liquid (Fig. 

Current understanding [31 [H] attributes the jerks to the mechanical buckling or wrinkling 
instability of any elastic sheet that is floating on a liquid and is under compression [27] . 
Above a threshold pressure Pc, the sheet distorts, at a cost of bending energy. The wrinkles 
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also produce a net upward and downward displacement of the liquid, thus increasing its 
gravitational energy. Taking account of these costs, one finds a threshold pressure given 
by Pc = '^y/B'pg, where B is the bending stiffness and g is the acceleration of gravity. The 
predicted buckling wavelength A is given by A = 2TT[B/{pg)Y/*. This buckling is unstable, 
leading to a release of the overpressure pc over a displacement Af comparable to A. Though 
Pc is too small to measure directly in our system, it can be estimated either using typical 
values of B or by inferring B from the observed jerk displacement. Bending stiffness for 
lipid monolayers like our system lie in the range B < 10^^^ J. [28j. This B value implies 
an upper bound for p^. Pc < Q x 10^^ mN/m. If instead we infer pc from the observed jerk 
displacements via pc = pgX^/{2n^) = pgA//{2TT^), we obtain values of 2 x 10"^ mN/m, 
i.e. 30 times smaller than our upper bound [29j. In what follows we will use the more 
conservative "upper-bound" value of pc, namely 6 x 10^^ mN/meter. 

To summarize the above discussion, the sheet is assumed to deform in the x-z plane 
while remaining uniform along the y axis. It is laterally compressed by a two-dimensional, 
uniform pressure p. (In the experimental system p is given by the actual pressure exerted at 
the boundaries minus the surface tension of the liquid.) Prior to instability the monolayer 
remains flat {h = 0) and responds to compression by slightly decreasing its actual length. 
When the pressure exceeds pc, the monolayer wrinkles or folds out of the x-y plane {h ^ 0). 
Since the pressure pc produces negligible elastic compression, the total length L is fixed, and 
the displacement A along the x axis is fully accounted for by the wrinkles or folds. The 
underlying liquid, having viscosity r] and mass density p, occupies the region z < h. 




FIG. 2. Schematic view of the system and its parameterization. 
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III. DYNAMICS 



A. Constant Pressure 



As noted above, the motion of the sheet depends strongly on the viscous entrainment 
of the fluid beneath it. In this section we determine the motion taking account of this 
entrainment. We begin with a simplified situation. Suppose we control the surface pressure 
and increase it above the instability pressure, pc- We further suppose that as the sheet 
begins to fold, the buckled region is decreasing the pressure by a constant amount p. At this 
moment the forces on the sheet are unbalanced and it begins to accelerate. In all but the 
small folding region the sheet is horizontal and its motion is a pure horizontal translation. 
The fiat sheet adjacent to the fold will translate laterally, thus creating a velocity profile 
v{z)Sc in the underlying liquid. The resulting viscous drag decreases with time, as the velocity 
profile penetrates deeper into the liquid and its gradient becomes less sharp. We shall assume 
that the sheet's velocity is always equal to that of the water immediately beneath it. That 
is, the fluid obeys a no-slip boundary condition. We examine such assumptions in Section 



IV below. We may neglect the mass of the sheet relative to the much greater mass of the 
entrained fluid. Thus the problem is simplified to that of finding the velocity of a half infinite 
fluid due to a rigid surface moving with a constant lateral pressure. 
The equation of motion for the fluid is: 

V = ud^zV, (1) 

where u is the kinematic viscosity given by the ratio of viscosity and density, u = rj/p, and a 
dot denotes derivative with respect to time. The constant stress at z = implies a boundary 
condition of the form 

dzv\z=o = (2) 
Accompanying it are the condition of vanishing flow far away from the interface, 

v{z ^ -oo,t) =0, (3) 

and initial conditions of stationarity, 

A{t = 0) = A{t = 0) = v{z, t = 0) = 0. (4) 
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Eqs Q and (g may be simplified by defining the Laplace transform v{uj) by v{u) = 
J,^v{t)e~-'dt: 

^ p 

uv{z,uj) = ud^zv{z,uj), d^v\^=o = — . (5) 

The solution in Laplace space is 



V[Z, u) 



and in the time domain, 



exp(yavV z), (6) 



v{z,t) = 2j9(z/t/(7rLV))'^'[e~°' + a(Erf(«) + 1)], a = z/{2V^t) (7) 

The displacement of the sheet is given by integrating the velocity on the surface, A = 
J^v{z = 0,t')dt', 

A = t^/^ (8) 

3v^ Lr] ^ ^ 



Since the accelerated mass is constantly increasing in time, the acceleration v decreases with 
time as This leads to the unusual "jerky" increase of the velocity v{t). We note that 

this predicted motion arises from basic hydrodynamics; it does not depend on the source of 
the pressure p. We note also that a fixed fraction ~ 0.64 of the input power pA goes into 
kinetic energy; the remainder goes into viscous dissipation. 



B. General Pressure 



Buckling relaxes stress at the tip of the fold; we would therefore expect the driving 
pressure to depend on the displacement. The equation for A{t) can readily be generalized 
to include such a p{A). We discuss it here for completeness. 

La A = p{A) - Lr]d, v\^^q, (9) 

where the first term on the right is the driving stress, the last term is the stress from the 
fluid, and a is the two-dimensional mass density of the sheet. For the fluid, Eq. ([T]) still 
applies, and the two are supplemented by a no-slip boundary condition at the interface, 

v{z = 0,t) = A, (10) 

a condition of vanishing flow far from the surface (Eq. [s]), and initial conditions of station- 
arity (Eq. 
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We proceed by eliminating v to obtain an equation for the sheet alone. The velocity 
Green function of the liquid is given by the boundary condition v{z = 0,t) = 6{t) and Eqs. 
( |3pl ). Explicitly, 



The general solution with the no-slip boundary condition of Eq. ( 10 ) is obtained simply 
by integration, v{z, t) = Jq G{t - t')A{t')dt' = Jq G{t')A{t - t')dt' . Inserting this expression 
in Eq. ^ and integrating by parts, produces an equation for the sheet alone: 

LaK=p[A)-^fdt''^^. (12) 



'TIV Jo 

As before, we may neglect the inertial term on the left hand side. For a constant pressure, 
Eq. (12) gives the same displacement as Eq.([7]). 



C. Application to observed jerks 

We may now compare the observed jerk motion with the motion expected from the folding 



forces. It is straightforward to solve Eq. (12) using the predicted (quadratic) p{A) [30]. The 
resulting motion proves to be much slower than the observed jerks. We may quantify this 
discrepancy in two ways. First we ask what constant pressure p would be required to give 
the observed displacements A in the observed time t. Then we determine the observed 
kinetic energy and compare it to the energy available from folding. 

To find the pressure required to produce the observed jerks, we substitute A = 2 microns 
and t = 0.12 sec into Eq. (12), to obtain p = 3.2 x 10~^ mN/m. This is some 50 times larger 



than the upper-bound folding pressure Pc = Q x 10~^ mN/m obtained above. 

The same discrepancy emerges if we compare the observed kinetic energy E^. with the 
work Wp done by the pressure difference p(A). Both are proportional to the width w of the 
jerking region. The pressure p(A) is always smaller than pc throughout the folding. Thus 
Wp/w < Pc A < 1.2 X 10"^^ N . Other things being equal, the kinetic energy Ek for a 
given average velocity A/t is larger if t is larger (since larger t implies a greater entrained 
mass). Thus we may obtain a lower bound on the kinetic energy by limiting t to the observed 
duration of 0.12 sec. For the moment we simply use the constant-pressure solution of Eq. ([t]) 
and use 



Ek/w = L -p J dz v{zy = 1.7 Ly/ui 



1 /A 



P > t 



(13) 
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We examine this estimate in Sec. IV For L = 0.05m, A = 2 microns and t = 0.12 sec, 
this gives Ek/w = 4 x 10~^^ N — some 32 times our upper bound of the work Wp supphed 
by folding. Thus even when one ignores the work that must go into viscous dissipation, the 
observed energy is far larger than what the folding energy ,Wp, can supply. (Taking account 
of the dissipated energy, one recovers the factor-50 discrepancy quoted above.) 



IV. DISCUSSION 

The arguments above indicate a worrisome discrepancy between the observed jerking 
motion and the mechanics of folding presumed to account for this motion. The discrepancy 
is serious; it survives even when we use conservative bounds in our estimates. In this section 
we survey possible ways to account for the discrepancy. First we review possible flaws in 
our description of the kinetic energy and the folding forces. Then we consider other forces 
that might account for the jerking motions. 

Our estimate of the kinetic energy was a simplified one, but it gives a proper lower 
bound for a given displacement A and time t. We will proceed by reinforcing several of 
the assumptions. Firstly, we assumed a constant pressure, though the actual unbalanced 
pressure increases with time. However other choices would have led to a higher kinetic 
energy. We consider the effect of replacing our constant-pressure estimate by allowing the 
pressure to increase with time. In order to achieve the required A(= J dt A) in the given 
time if: with a time-increasing pressure, we will necessarily reduce A at early times and 
increase it at late times. However, any shift of A from earlier to later times has the effect of 
increasing the kinetic energy. To see this, we consider a small decrease of A at time t< over 
a brief interval At. To maintain a fixed total displacement A, we make an equal addition 
to A at a later time t>. Any shift of A from earlier to later times can be accomplished by 
repeating this process. This perturbation of A creates a corresponding perturbation of the 
fluid velocity at the final time t: we denote it by 6v{z,t). We may then express the final 
kinetic energy Ek{t) in terms of this 6v and the initial profile vo{z,t) using the integral of 
Eq.[7} 

Ek{t)/w = L^pJdz [vo{z,t) + 6v{z,t)]^ 

= Eok{t)/w + LpJdz vo{z, t) 6v{z, t) + 0{Sv^) (14) 
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The second perturbing term is necessarily positive provided the (positive) vq{z) profile is 



monotonic. To see this we express 6v in terms of the Green function G of Eq. (11): 



5v{z, t) = 5 At [-G{z, t - t<) + G{z, t - t>)]. (15) 

We may express any monotonic Vq as a sum of positive step functions extending from 
to some Z. For a given step function, the contribution to is given by 6v. Thus it 
suffices to show that this integral is positive. In tecarms of the G functions, this means 



G{z, t — u) is an increasing function of u. This may be verified explicitly using Eq. ( 11 ). 
Thus a small shift in the pressure profile from a constant one to an increasing one with 
the same A and t only increases the energy E^. If further small shifts are added, the same 
reasoning implies that the E^ again increases, provided the starting Vq{z) remains monotonic 



in z. We conclude that the constant-pressure Ek of Eq. (13) under-estimates Ek- 

A second assumption that affects our estimate of the kinetic energy is the no-slip boundary 
condition between the sliding monolayer and the fluid beneath. Having a slip would cause 
less drag of the fluid and thus less kinetic energy. But in our case of hydrophilic heads facing 
the water, there is no justiflcation for a signiflcant slip. 

A further effect that can potentially reduce the kinetic energy is the possibility that the 
motion is restricted in area, so that less fluid is entrained. As noted above, this restriction 
can occur if the surface layer is compressible. The expected compressibility is such that there 
could be a noticeable effect on the scale of the entire 15 cm sample. To account for this, 
we assumed that the motion is restricted to a range of only 5 cm. On this scale we argued 
that any departure from rigid sliding of the sheet would be negligible. We conclude that 
our estimate of the kinetic energy is indeed a lower limit, given the observed displacements 
A and times t. Thus the factor-32 discrepancy between the kinetic energy and the folding 
pressure remains. 

Another potential way to resolve this discrepancy is that we have txndenestimated the 



folding pressure p. As was noted in section III, the creation of a fold reduces the pressure 



It cannot possibly be larger than the instability pressure, Pc = '^^/Bpg that we assumed 
above. Would additional folds increase the available folding energy? No. The energy must 
be smaller than pc times the observed displacement, as we assumed. Could it be that 
> 6 X 10^^ mN/m? Only if we underestimated the bending rigidity, B. The local 
structure of the monolayer might imply a bending rigidity greater than is usually observed 
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due to texture in the sheet [3T], but how much larger? Equating the folding and kinetic 
energies requires Pc of at least 2 x lO^^N/m i.e. B > 10^^^ J — a suspiciously larger value. 

The kinetic energy could not be smaller than what we have assumed, and folding energy 
is by no means larger, so the inconsistency stands. It can only be resolved if other forces 
are considered. Many ambient forces in the experiment could easily produce pressures 
larger than pc- Examples are the measured surface pressure, the residual surface tension, or 
ambient vibration. However, these ambient effects are not adequate to explain the observed 
folding instability. In order to cause an instability, the unbalanced pressure must increase as 
the displacement increases. The ambient pressures listed above do not have this property. 
Another collapse mode, e.g. brittle fracture, would indicate that a different mechanism is 
taking place, but the folding was observed to be a reversible process as opposed to fracture. 

Another potential destabilizing force is adhesion. There is an energy gain when two 
hydrocarbon surfaces are in mutual contact. The adhesion energy per unit area is of order 
2 X 10~^J/m^, more than enough to resolve the entire discrepancy. Adhesion energy was 
invoked to account for other, fluid-like collapsed structures [lOj, but was not considered in the 
accepted explanation for the straight abrupt folding [11]. From the arguments given above 
such an effect appears to be necessary, without which jerking would not be observed. We 
will note that adhesion of the hydrocarbon chains requires a large, systematic, deformation 
of the sheet prior to jerking. In that case there is no longer an obvious relation between the 
length of scale for wrinkles and folds and the length scale for jerk displacement. 

The force responsible for jerks in the monolayers of Ref. [13] may well be important 
in a broader context. The wrinkle-fold transition has been implicated in a broader class 
of nanoscale systems: nanoparticle trilayers [32] and single-component lipid monolayers [3|. 
The various energies in these systems are different, and the motion is also somewhat different 
from that of Ref. [13] . Nevertheless, the discrepancy shown above may also apply to these 
other systems. It is also possible that the force responsible for jerks affects structure as 
well as dynamics. That is, the magnitude of a fold may be governed by other characteristic 
lengths than the wrinkle wavelength A. Still, the incipient instability might be due to the 
wrinkle-fold mechanism while the rapid subsequent motion is controlled by adhesion or some 
other force. 

This work highlights the distinctive dynamics of thin solid sheets on fluids in general. 



Equation (12) apphes whenever such a sheet is accelerated by an unbalanced force, whether 



13 



folding occurs or not. The equations imply a distinctive form of acceleration and a distinctive 
partitioning of the supplied work into kinetic energy and dissipation. 

V. CONCLUSION 

The wrinkle-to-fold mechanism [Tl] for thin film buckling has allowed a new avenue 
for understanding buckling phenomena in a range of nanoscale systems. Its successes in 
explaining structural aspects of folding have led us to apply it to the well-studied dynamics 
of monolayer jerks. For these dynamic phenomena, our study indicates that other forces are 
at play. Supporting this conclusion are experiments on lipid-coated micro-bubbles, where 
it is clear that gravity, an important ingredient of the folding mechanism, plays no role 
fl2\ [15] . These other forces must be much stronger than those previously considered. The 
necessity of such forces underscores the remarkable nature of monolayer jerks. It also raises 
the importance of understanding these jerks. 
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